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ABSTRACT 

We report the first detection of the Integrated Sachs- Wolfe effect (ISW) in wavelet 
space, at scales in the sky around 0^7° with a significance « 3.3cx, by cross-correlating 
the Wilkinson Microwave Anisotropy Probe (WMAP) first-year data and the NRAO 
VLA Sky Survey (NVSS). In addition, we present a detailed comparison among the 
capabilities of three different techniques for two different objectives: to detect the 
ISW and to put constraints in the nature of the dark energy. The three studied tech- 
niques are: the cross-angular power spectrum (CAPS, harmonic space), the correlation 
function (CCF, real space) and the covariance of the Spherical Mexican Hat Wavelet 
(SMHW) coefficients (CSMHW, wavelet space). We prove that the CSMHW is ex- 
pected to provide a higher detection (in terms of the signal-to-noise ratio) of the ISW 
effect for a certain scale. On the other hand, the detection achieved by the CAPS is 
the lowest one, being the signal-to-noise ratio dispersed among a wide multipole range. 
The CCF provides an intermediate detection level. This prediction has been corrobo- 
rated by the analysis of the data. The SMHW analysis shows that the cross-correlation 
signal is caused neither by systematic effects nor foreground contamination. However, 
by taking into account the information encoded in all the multipoles/scales/angles, 
the CAPS provides slightly better constraints than the SMHW in the cosmological 
parameters that define the nature of the dark energy. The limits provided by the CCF 
are wider than for the other two methods, although the three of them give similar 
confidence levels. Two different cases have been studied: 1) a flat ACDM universe and 
2) a flat universe with an equation of state parameter that, although does not change 
with time, could take values different from -1. In the first case, the CAPS provides (for 
a bias value of b = 1.6) = 0.73±„ " (at la CL). Moreover, the CAPS rejects the 
range fl\ < 0.1 at w 3.5er, which is the highest detection of the dark energy reported 
up to date. In the second case, the CAPS gives £Ide = 0.70^q' 2 o an d w = —Q.75i_fc\y 
(at la CL). This is the first estimation of the equation of state of the dark energy made 
through the cross-correlation of the CMB and the nearby galaxy density distribution. 
It also provides an independent estimation from the one made by the WMAP team 
using CMB and LSS. 

Key words: cosmic microwave background, cosmology: observations 



1 INTRODUCTION 

During the last decade, the study of the anisotropics of the 
cosmic microwave background (CMB) is being (probably) 
the most important source for understanding the universe. 
Since the first all-sky detection reported by COBE-DMR in 
1992 (Smoot et al. 1992), many ground and balloon based 



experiments have been measuring the CMB anisotropies, 
drawing (together with other astrophysical data like those 
provided by supernova) a flat universe whose evolution is 
dominated by the dark energy. This general picture has 
been recently confirmed by the first-year results of the 
NASA WMAP satellite (Bennett et al. 2003a, Spergel et al 
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2003), providing a strong confirmation of the fiducial ACDM 
model. 

At this point, one of the most interesting questions is 
the confirmation of this fiducial model by other indepen- 
dent analysis. In that sense, the Integrated Sachs- Wolfe ef- 
fect (ISW, Sachs & Wolfe, 1967) -that is produced by the 
time variation of the gravitational potential, in the linear 
regime- plays a crucial role since it provides either a direct 
indication of the presence of dark energy in the case of a 
flat universe or else the existence of spatial curvature (Pee- 
bles & Ratra 2003). Since WMAP has established strong 
constraints in the flatness of the universe (Bennett et al. 
2003a), the detection of the ISW will imply the detection of 
the dark energy. 

Many different cosmological models can be included 
within the dark energy concept: the standard inflationary 
model dominated by a cosmological constant (with an equa- 
tion of state parameter w = — 1; p = wp) and alternative 
models dominated by a dark energy whose energy density is 
spatially inhomogeneous with negative pressure and evolv- 
ing with time, like topological defects, quintessence models 
and phantom models (see for instance Melchiorri 2004 and 
references therein). However, there is not yet a convincing 
explanation for the origin and nature of the dark energy 
within the practical physics framework. Obviously, any in- 
formation extracted from cosmological data, that could help 
to characterise it, is very valuable (for phenomenological as- 
pects of dark energy see for instance Corasaniti 2004) . 

Crittenden & Turok (1996) firstly suggested that the 
ISW could be detected by cross-correlating the CMB with 
the nearby galaxy density distribution. The first attempt to 
detect the cross-correlation of the CMB and the galaxy den- 
sity distribution was done by Boughn & Crittenden (2002) 
using the COBE-DMR map and the NRAO VLA Sky Survey 
(NVSS, Condon et al. 1998) data. In that work, no cross- 
correlation was found, concluding that a future experiment 
with better sensitivity and resolution was required in or- 
der to detect it. The first-year WMAP data has provided 
a unique tool to perform such correlation. Several groups 
have reported ISW detection at different significant levels. 
Boughn & Crittenden (2004) have correlated the WMAP 
data with two different tracers of the nearby universe: the 
hard X-ray data provided by the HEAO-1 satellite (Boldt 
1987) and the NVSS data. They find a statistical signifi- 
cance of 1.8 - 2.8 a at scales lower than 3°. Nolta et al. 
(2004) have done an independent analysis of the WMAP 
and NVSS cross-correlation obtaining > at 95%. Fos- 
alba & Gaztanaga (2003) find a cross-correlation between 
WMAP and the APM Galaxy survey (Maddox et al. 1990) 
at 98.8% significance level at scales w 4° - 10°. Fosalba et al. 
(2003) cross-correlated WMAP with the Sloan Digital Sky 
Survey (SDSS, Abazajian et al. 2004) finding a positive ISW 
signal at the 3<r level and determining = 0.69 — 0.86 (2a 
CL). These two data sets were also analysed by Scranton 
et al. (2003)*. Afshordi et al. (2004) have cross-correlated 
the galaxies from the Near-IR Two Micron All Sky Survey 



(2MASS, Skrutskie et al. 1997) finding an ISW signal at 
2.5a level. 

All the previous works have performed the cross- 
correlation of the WMAP data and the nearby universe trac- 
ers using the cross-angular power spectrum (CAPS, e.g. Af- 
shordi et al. 2004) and the correlation function (CCF, e.g. 
Boughn & Crittenden, 2002, 2004, Nolta et al. 2004). In the 
present paper, we propose an alternative to the harmonic 
and real spaces to perform the correlation of the WMAP 
and NVSS data sets: the wavelet space. The motivation for 
using wavelets is clear since the ISW signal provided by the 
CMB-galaxies cross-correlation takes place at scales between 
2° — 10° (Afshordi 2004). As it is well known, wavelets are 
very suitable tools for detecting signals with a characteristic 
scale: by filtering the data at a given scale, those structures 
having that scale are amplified, allowing for a most opti- 
mal detection. Therefore, the wavelet analysis of the cross- 
correlation will concentrate almost all the signal of the ISW 
detection in a small scale range, allowing for a fast and direct 
measure of its signal-to-noise. This is not possible neither for 
the CAPS nor the CCF, where all the information spread 
out among all the multipoles/angles must be considered to 
estimate the total signal-to-noise. This effect is even more 
important for the CAPS than for the CCF. 

Since we expect that the common structures in CMB 
and galaxy surveys have roundish rather than elongated 
shapes, we propose the Spherical Mexican Hat Wavelet 
(SMHW) to study the cross-correlation. A natural statis- 
tic of the cross-correlation in wavelet space is the covariance 
of the SMHW coefficients (CSMHW) at each scale. Wavelets 
have been extensively used in many astrophysical and cos- 
mological analyses. In the CMB field, they have been used to 
study the non-Gaussianity of the CMB (Pando et al. 1998, 
Hobson et al. 1999, Aghanim & Forni, 1999, Tenorio et al. 
1999, Barreiro et al. 2000, Mukherjee et al., 2000, Barreiro 
& Hobson 2001, Cayon et al. 2001, 2003, Matmez-Gonzalez 
et al. 2002, Vielva et al. 2004, Starck et al. 2004, Cruz et 
al. 2005, McEwen et al. 2005), the component separation 
problem (Cayon et al. 2000, Vielva et al. 2001a, b, 2003) 
and denoising (Sanz et al. 1999a, b). This work is the first 
application to the cross-correlation of CMB and other data 
sets. 

In addition to the pure ISW detection, the CAPS, the 
CCF and the CSMHW can be exploited, by comparing the 
ISW signal with the one provided by different cosmological 
models, to put constraints in the nature of the dark energy: 
not only in the quantity of the dark energy (Qde), but also 
in the equation of state parameter (w). 

The paper is organized as follows. In Section 2 we de- 
scribe the analysed data sets (WMAP and NVSS). In Sec- 
tion 3 the three cross-correlation estimators (CAPS, CCF 
and CSMHW) are presented. The expected signal-to-noise 
ratio for the ISW detection for the CAPS, the CCF and the 
CSMHW is calculated in Section 4. The results are given in 
Section 5 and, finally, in Section 6 are the conclusions. 



* See Afshordi et al. (2004) for a critical discussion of the analysis 
done in Fosalba & Gaztanaga (2003), Fosalba et al. (2003) and 
Scranton et al. (2003) 



2 THE WMAP AND NVSS DATA SETS 

The two data sets that have been used in order to perform 
the CMB-nearby universe cross-correlation are the Wilkin- 
son Microwave Anisotropy Probe (WMAP, Bennett et al. 
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Figure 1. Analysed WMAP and NVSS data after the application of the joint mask and the subtraction of the residual monopole and 
dipole. The maps are represented in the HEALPix scheme, with a resolution parameter N si d e = 64 pixel size fs 55 arcmin). 



2003a and references therein) first-year data and the NRAO 
VLA Sky Survey (NVSS, Condon et al. 1998). 

2.1 WMAP data 

The WMAP radiometers observe at five frequencies: 22.8, 
33.0, 40.7, 60.8 and 93.5 GHz, having 1, 1, 2, 2 and 4 
receivers per frequency band, respectively. All the papers, 
data and products generated by the WMAP team can be 
found in the Legacy Archive for Microwave Background 
Data Analysis (LAMBDA) Web sitet. The WMAP maps 
are presented in the HEALPix scheme (Gorski et al. 2005) 
at the N S id e = 512 resolution parameter. The WMAP team 
and other groups have proposed different CMB maps ob- 
tained from the WMAP data. In this work, we have used 
the map proposed by the WMAP team and already used 
by other groups (Komatsu et al. 2003, Vielva et al. 2004, 
Eriksen et al. 2004, Mukherjee & Wang 2003, Hansen et al. 
2004, Cruz et al. 2005) for the study of the non-Gaussianity 
and the isotropy of the CMB. This map is generated (see 
Bennett et al. 2003b for details and Vielva et al. 2004 for 
a summarized description) as the noise weighted combina- 
tion of all the maps produced by the receivers in which the 
CMB is the dominant signal (40.7, 60.8 and 93.5 GHz), 
after subtraction of the foreground emission and applica- 
tion of the so-called "KpO" Galactic mask (defined by the 
WMAP team and where the brightest point sources are also 
masked). Whereas for the non-Gaussianity studies the res- 
olution N si de = 256 was commonly chosen, in the present 
work we have degraded the combined, corrected and masked 
map down to N s ide = 64 (pixel size « 55 arcmin). The rea- 
son is that, as it was pointed out by Afshordi (2004), almost 
all the signal of the ISW is expected to be generated by 
structures with a scale larger than 2° . Hence, a WMAP res- 
olution of around 1° is enough. 

2.2 NVSS data 

The NVSS catalogue covers around 80% of the sky and has 
flux and polarization measurements for almost 2 million of 
point sources with a minimum flux ~ 2.5 mjy at 1.4 GHz. 

t Available at http://cmbdata.gsfc.nasa.gov 



This catalogue has been already used for performing the cor- 
relation with the WMAP data (Boughn & Crittenden 2002, 
2004 and Nolta et al. 2004). In this work, we have repre- 
sented the point source catalogue in the HEALPix scheme, 
also at the N a ide = 64 resolution. Only the sources above 2.5 
mjy have been used, which represents the 50% completeness 
(Condon et al. 1998). As it was pointed out by Boughn & 
Crittenden (2002) and Nolta et al. (2004), the mean density 
of point sources varies as a function of the declination. This 
systematic effect was corrected, by imposing that the mean 
of the galaxy density at each iso-latitude band is the same. 
The iso-latitude bands are defined taking into account the 
change in the RMS noise levels in the NVSS survey (see 
Figure 10 of Condon et al. 1998). Alternative strategies (as 
the one proposed by Nolta et al. 2004) were also consid- 
ered, proving that results are not sensitive to the particular 
correction procedure. 

As it was said before, the NVSS catalogue covers around 
80% of the sky: for equatorial declination lower than —50° 
there are not observations and, within the range —37° > 
S > —50° the coverage is not good enough. Hence, we only 
consider sources with an equatorial declination 5 > —37°. 
With all these constraints, we have a galaxy distribution 
map of ~ 1600000 radio sources with an average number of 

40.4 counts per pixel. 

In Figure 1 we have plotted the two maps to be anal- 
ysed: WMAP (left) and NVSS (right). Both are in Galactic 
coordinates and the joint mask (KpO + S < —37°) has been 
applied. The residual monopole and dipole outside the mask 
have been removed. 

2.3 Simulations 

We have also performed realistic simulations in order to 
carry out the analysis. 1000 Gaussian simulations of the 
WMAP data have been done, following the concordance cos- 
mological model given by the Table 1 of Spergel et al. 2003 
(fi A = 0.71, On = 0.29, Q b = 0.047, H = 72, r = 0.166, 
n = 0.99) and using the CMBFAST code (Seljak & Zaldar- 
riaga 1996). For each realization, we have simulated all the 
WMAP data measured by the receivers at 40.7, 60.8 and 

93.5 GHz, they have been convolved with the real beams 
provided in the LAMBDA Web site, the anisotropic WMAP 
noise was added, the maps were combined using a noise- 
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weighted average, the combined map was degraded to the 
N S ide = 64 resolution and the joint mask was applied. We 
have cross-correlated the 1000 CMB simulations with the 
NVSS data, in order to evaluate the significance level of the 
cross-correlation obtained from WMAP and NVSS. This is 
enough to quantify the covariance matrix associated to ran- 
dom cross-correlations and we have checked that it is almost 
independent of the cosmological model used to simulate the 
CMB. 



3 THE CROSS-CORRELATION ESTIMATORS: 
CAPS, CCF AND CSMHW 

As it is well known (e.g. Peebles & Ratra 2003) for a flat 
universe where the dynamics are dominated by the dark 
energy, we expect a positive correlation between the CMB 
and the galaxy distribution of the nearby universe (z < 1, 
see for instance Afshordi 2004). 

In this paper we apply three different techniques to 
study such correlation and to compare the performance 
of each of them for the detection of the ISW. The three 
studied techniques are the cross-angular power spectrum 
(CAPS), the correlation function (CCF) and the covari- 
ance of the SMHW coefficients (CSMHW), covering the har- 
monic, real and wavelet spaces. The CAPS has been already 
used, for instance, in Afshordi et al. (2004) to estimate the 
Sunyaev-Zeldovich, the point sources and the ISW signals 
by cross- correlating the 2MASS infrared source catalogue 
and WMAP. The CCF has been more extensively used, for 
instance in Boughn & Crittenden (2004) and Nolta et al. 
(2004). We proposed a new technique based in wavelet space, 
the CSMHW. 

3.1 The cross-angular power spectrum 

The observed CAPS of WMAP and NVSS is given by: 

CewN = (ti m n^ m ) , (1) 

where te m and ne m are the harmonic coefficients of the 
WMAP and the NVSS maps (respectively) and * denotes 
complex conjugation. 

This CAPS can be compared with the theoretical one 
expected from a given cosmological model: 

C\ heo wN = Pe 2 bwe oni M til i C e i the ° WN , (2) 

where pi is the pixel window function for the corresponding 
pixelization, bwe and bne are the beam window function for 
WMAP and NVSS*, respectively, M e ji is a coupling matrix 
that accounts for the joint mask (see Hivon et al. 2002) and 
C/ heo wN is the pure CAPS of WMAP and NVSS. Since 
N a ide = 64 resolution maps are analysed, only multipoles 
from £ = 2 to £ — 191 are considered. As it can be seen in, 
e.g., Nolta et al. (2004), C e theo wN is given by: 

C e the ° WN = \2M m H 2 J ^A s \k)Ft W (k)F e N (k), (3) 

where Q. m is the matter density, H is the Hubble constant, 

t For the NVSS galaxy density map there is no beaming, there- 
fore b Ne = 1 V I 



As (k) — k Pg(k)/2ir is the logarithmic matter power spec- 
trum (being Ps(k) the matter power spectrum) and Fe w (k) 
and Fe N (k) are filter functions for the CMB and the galaxy 
density distribution, respectively, given by: 

Fi W (k) = J dzf z j e (kr,(z)) (4) 

F e N (k) = b J dz^-D(z)j e (kn(z)) (5) 

In the case of Fe w (k), the integral is defined from z = to z 
at recombination epoch. The integration range for the filter 
function Fe N (k) is defined, in practice, by the source redshift 
distribution function The function D(z) is the linear 
growth factor for the matter distribution (calculated from 
CMBFAST by computing the transfer function for different 
redshifts), g = (1 + z)D(z) is the linear growth suppression 
factor, je(kr}(z)) is the spherical Bessel function and 7)(z) is 
the conformal look-back time. The bias factor b is assumed 
to be redshift independent. 

For the evolution of D(z) we consider the standard 
model dominated by a cosmological constant as well as al- 
ternative models dominated by a dark energy whose energy 
density is spatially inhomogeneous with negative pressure 
and evolving with time, p = po(l + z) 3 ^ 1+ ™\ where w is 
the equation of state parameter defined as the ratio of pres- 
sure to density. Whereas the standard inflationary model 
assumes w = — 1 and a homogeneous field (the cosmological 
constant), in general dark energy models are characterised 
by values of w < 0. For instance, topological defects can be 
phenomenologically represented by an equation of state pa- 
rameter — 1/3 > w > —2/3 (see e.g. Friedland et al. 2003); 
quintessence models (Wetterich 1988, Caldwell et al. 1998) 
imply —1 < w < and phantom models have equation of 
state parameters w < — 1 (these last models, however, vio- 
late the null dominant energy condition; see Carroll et al. 
2003 for a detailed discussion). 

Although a convincing explanation for the origin and 
nature of the dark energy within the framework of particle 
physics is lacking, the models considered in the literature 
produce in general w varying with time. However, for most 
dark energy models, the equation of state changes slowly 
with time and a standard approximation assumes that (at 
least during a given epoch) w can be considered as an effec- 
tive (and constant) equation of state parameter (Wang et 
al. 2000). Therefore, in this paper, we consider w constant 
as an useful approach to extract fundamental properties of 
the dark energy. 

3.2 The correlation function 

The observed CCF is given by: 

CCF(8) = 1 g T(n)N(n'), (6) 

where T(n) and N(n) are the WMAP temperature and 
NVSS galaxy density (respectively) at position n. The sum 
is extended over all the pixels that are allowed by the joint 
mask and inside a disk described by (9 — A6/2, 6 + A9/2) 
where nn" = cos(6) and A6 = 2°. As in Nolta et al. (2004) 
the CCF is evaluated at 10 angles 9 = 1, 2, 3, 5, 7, 9, 10, 
13, 15, 17 and 19 degrees. 
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It can also be compared with the theoretical one ex- 
pected from a given cosmological model: 



CCF theo (9) = (ST(n)6N(n) 



(7) 



2 , + 1 Pe(nn) p 2 b W e b N i C/ 



where Pe(nh") are the Legendre polynomials and Ce the ° WN 
is given by Equation 3. 

3.3 The covariance of the wavelet coefficients 

We define the CSMHW at a given scale R as: 

Cov WN (R) = ^-^u }T (R,p)iv N (R,p), (8) 

V 

where lu t (R, p) and lj n (R, p) are the SMHW coefficients 
of WMAP and NVSS (respectively) at position p. We pro- 
pose to study the cross-correlation of the WMAP and NVSS 
Spherical Mexican Hat Wavelet (SMHW) coefficient maps 
at different scales (i?i = 50.0, R 2 = 75.0, R 3 = 100.0, 
Ri = 150.0, R s = 200.0, R e > = 250.0, R 7 = 300.0, R» = 
400.0, R 9 = 500.0, R 10 = 600.0, Ru = 750.0, R 12 = 900.0, 
R13 = 1050.0 arcmin) in the angular range where the ISW 
signal is expected to be more important (Afshordi 2004) for 
CMB-galaxy cross-correlations. 

The SMHW coefficients are obtained by convolving the 
map with the SMHW at a given scale: 



oj x (R,p) = J dn'x(p + p')^ s (e',R), 

where <5>s(0',R) is the SMHW: 



(9) 



Vy ' ' ^N{R)\ 
being N(R) a normalization constant: 

R 2 R^ 1/2 



-y 2 /2R 2 



N(R)^R(l + ^ + ^) 



(10) 



(11) 



The distance on the tangent plane is given by y that is re- 
lated to the latitude angle (6) through: 



y = 2 tan -. 

y 2 



(12) 



The SMHW is obtained from the Euclidean Mexican Hat 
Wavelet (MHW) following the stereographic projection sug- 
gested by Antoine & Vanderheynst (1998). The reader is re- 
ferred to the work of Martinez- Gonzalez et al. (2002) where 
this projection is described for the SMHW, as well as its 
properties. 

In Equation 8, the sum J2 p ls extended over all the pix- 
els that are not masked (Nr). As it was already discussed 
in Vielva et al. (2004), since we are convolving a map with a 
large mask (the joint mask covers 43% of the sky) with the 
SMHW at a given scale R, we are introducing into the cross- 
correlation estimator Cov WN (R) a large number of pixels - 
those close to the border of the mask- highly affected by the 
zero value of the mask, which implies a loss of efficiency. For 
that reason, at a given scale R, we exclude from the calcula- 
tion of Cov WN (R) all the pixels with a strong contamination 
from the joint mask. Whereas in Vielva et al. (2004) those 



pixels that were closer than 2.57? to any one of the pixels of 
the joint mask where excluded, for the present analysis we 
exclude only those pixels closer than 1.0R to the joint mask. 
There are two reasons for that: first, the mask used in the 
present work is much larger than the one used in Vielva et 
al. (2004) and second, the number of pixels is lower, since 
the HEALPix resolution is 3 times lower. These two facts 
make the available number of data much smaller (for all the 
scales) than the ones in Vielva et al. (2004) and, therefore, 
we must be less restrictive in order to construct the 13 exclu- 
sion masks. However, we have checked that the ISW is also 
detected with the more restrictive exclusion masks proposed 
by Vielva et al. (2004). 

The curve defined by Equation 8 can be easily compared 
with the theoretical CSMHW between the CMB and galaxy 
density distribution: the Cov WN (R) is nothing more than 
the mean value of a map M(R) — lu t (R)lo n (R), where the 
reference to the position p has been removed for simplicity. 
It is straightforward to show that the theoretical wavelet 
cross-correlation can be written as: 



Cov t 



(13) 



where Ce M (R) is the angular power spectrum of the map 
M(R) and it is given by: 

Ce M (R) = Pe 2 s 2 e (R) b wi b Nl C e theo WN , (14) 

where, as before, pe is the pixel window function for the cor- 
responding pixelization, b W e and b Ne are the beam window 
function for WMAP and NVSS, C t theo WN is the pure CAPS 
of WMAP and NVSS and se(R) are the spherical harmonic 
coefficients of the SMHW at scale R (there is no dependence 
with m, since the SMHW is isotropic). 

Summarising, the comparison of the experimental curves 
given by Equations 1, 6 and 8 with the theoretical prediction 
given by Equation 2, 7 and 13 can be done, by computing 
the filter functions and the pure CAPS for different cosmo- 
logical models. 



4 COMPARISON OF DIFFERENT 

TECHNIQUES: MOTIVATION FOR USING 
WAVELETS 

The motivation for using wavelets for the detection of the 
ISW signal through the CMB-galaxies cross-correlation is 
given by the fact that such correlation has his maximum 
contribution in the scale range between 2° and 10° (see Af- 
shordi 2004). As it is well known, wavelets are ideal tools 
for detecting signals (embedded in a given background) with 
a characteristic scale. The convolutions of the data with a 
wavelet with a certain scale, amplifies those structures with 
such scale, allowing for a most optimal detection. 

In the present Section, we compare the expected signal- 
to-noise ratio of the cross-correlated signal between CMB 
fluctuations and the radio sources density distribution as a 
function of the angular scale for the estimators described in 
the previous section (see Figure 2): the cross-angular spec- 
trum Ce WN (CASP, Eq. 3); the correlation function (CCF, 
Eq. 7) and the covariance of wavelet coefficients Covwn (R) 
(CSMHW, Eq. 13). This is done for an ideal free-noise ex- 
periment with full-sky coverage and a resolution N Bide = 64 
(and assuming no beaming, i.e. bwe = bNe = !)■ 
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Figure 2. The expected signal-to-noise ratio of the CMB-radio 
sources cross-correlation according to the three studied tech- 
niques: the CSMHW (solid curve), the CAPS (dotted-dashed 
curve) and the CCF (dotted curve). The cosmological parame- 
ters are chosen according to the concordance model. 



Let us recall two important points. First, the values of 
the signal-to-noise ratio curves for both, the CSMHW coef- 
ficients and the CCF have a very high scale-to-scale correla- 
tion, whereas there is no correlation in the case of the CAPS. 
In principle, by taking into account all the scales and the full 
correlations, one should expect to get comparable detection 
levels in the three cases. However, we remark that, for a cer- 
tain scale, the detection of the ISW given by the CSMHW 
is proven to be the highest. In other words, the whole in- 
formation of the ISW has been extremely concentrated and 
amplified in a very narrow range of angular scales. This is 
undoubtedly very useful in a practical situation, like the one 
addressed in this paper. 

Second, in the case of a noisy experiment, the signal- 
to-noise ratio should be lower than the one obtained under 
our ideal conditions. This will be much more important for 
an estimator as the CCF, since it has its maximum amplifi- 
cation at the smallest scale. On the contrary, the CSMHW 
and the CAPS are almost unaltered in this case, since the 
maximum amplification is reached at intermediate angular 
scales, where the noise is almost negligible. The same argu- 
ment is also applied if beaming is included. 



5 RESULTS 



The signal-to-noise level can be estimated as the ratio 
between the cross-correlation estimator and its dispersion. 
As it is straightforward to show, the dispersion of the CAPS 
for two Gaussian random fields is given by: 



Ctheo^ i /~i theo 1-1 theo 



21 + 1 



(15) 



where Ce the ° w and Ct theo N are the angular power spectra 
for the CMB and the radio sources density distribution, re- 
spectively. 

Using Equation 13 and 15, we find that the dispersion 
of the correlation of the CSMHW is given by: 

AC*v the ° WN (R) = 



cw 



Gov 1 



E 2 4t^Pz A s t *(R)(C/ h '° 2 WN + C^\C t ^° w ) (16) 



16tt 2 



Equivalently, the dispersion of the CCF is given by: 



A.CCF eo (6) = 



(cCF th£ °{6) 2X j - (cCF th£ °{6)y 



E 



^^Pi 4 Pi 2 (nn>)(C t th ™ 2 WN + C t th »° N C t th °° w ) (17) 



Figure 2 shows that the CSMHW gives the highest 
signal-to-noise ratio (~ 3.5), at scales of « 4°. As expected, 
the maximum significance for the CCF (w 3) is reached at 
the smallest scale. The signal-to-noise ratio for the CAPS is 
significantly lower than for the other ones and always below 
1, with a shape nearly flat (ss 0.6). 



The results obtained from the cross-correlation of WMAP 
and NVSS are presented in this Section. In Subsection 5.1 we 
show the ISW detection at « 3.3cr provided by the CSMHW. 
The curves obtained for the CCF and the CAPS are also 
presented. In Subsection 5.2 we compare the obtained cross- 
correlation with different theoretical curves provided by dif- 
ferent ACDM cosmological models, in order to constrain 
the value of 57a- The performance of the three techniques 
(CAPS, CCF and CSMHW) is compared, showing that all of 
them provide compatible limits at la CL, being the CAPS 
and the CSMHW better than the CCF in order to reject 
a zero dark energy value. Both of then discard values of 
SI a < 0.1 with a significance > 3a. In particular, the CAPS 
detect the presence of dark energy at 3.5<r, which is the high- 
est detection reported up to date by using the ISW. Finally, 
in Subsection 5.3 we extend the analysis in order to con- 
strain the value of dark energy density and the equation of 
state parameter, finding again compatible values at la for 
the three statistics. 



5.1 The Detection of the ISW 

We have computed the SMHW coefficients of WMAP and 
NVSS maps at different scales, from 50 arcmin to 17.5° 
(which implies an approximate size in the sky from 1.5° 
to 35°) and we have studied their cross-correlation. Accep- 
tance intervals for the correlated signal at certain signif- 
icance levels a (32%, 5% and 1%) have been established 
at each scale using simulations. These acceptance intervals 
contain a probability 1 — a and the remaining probability 
is the same above and below the interval (i.e. a/2). The 
acceptance intervals have been determined by studying the 
CSMHW Cov WN {R) at each scale independently and calcu- 
lating the corresponding percentiles. 

In Figure 3 we present the ISW detection. Scales R$, 
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Figure 3. Cross-correlation in wavelet space of the WMAP and NVSS data (stars). Acceptance intervals for the 32% (inner interval), 
5% (middle interval) and 1% (outer interval) significance levels are also plotted. The mean value of the cross-correlation for the 1000 
simulations is also given (solid line). The plot is magnified for the lowest scales in the right panel. 



x 10" 3 
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2 4 6 8 10 12 14 16 18 10' 1Q 2 

6 (degrees) multipole (I) 

Figure 4. Cross-correlation in real (left panel) and harmonic (right panel) of the WMAP and NVSS data (stars). Acceptance intervals 
for the 32% (inner interval), 5% (middle interval) and 1% (outer interval) significance levels are also plotted. The mean value of the 
cross-correlation for the 1000 simulations is also given (solid line). 



f?6 and i?7 are outside the acceptance interval at 1% sig- 
nificance level (in particular, by increasing the number of 
simulations up to 10000, the scales Rs and R§ provide a 
detection at ~ 3.3<r), representing a clear detection of the 
cross-correlation of the WMAP and NVSS data. These scales 
correspond to a size in the sky around 6° - 8°, which is 
in very good agreement with the theoretical detection pre- 



dicted in Section 4 and also with Afshordi (2004), that pre- 
dicts an optimal scale for the ISW between 2° - 10°. 

For comparison, the detection curves for the CCF (left 
panel) and the CAPS (right panel) are given in Figure 4. 
The CCF provides a detection at ~ 95%, in good agree- 
ment with the curves obtained by Boughn & Crittenden 
(2004) and Nolta et al. (2004). For a better visualization, 
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Figure 5. The left panel shows the CMB-nearby galaxy density correlation along the WMAP frequency range (W, V, Q, K - 2.64Ka 
and combined WMAP maps). The right panel shows the covariance between the noisy W-V-Q map and NVSS. The acceptance intervals 
are the same as the ones in Figure 3 since they are dominated by the cosmic variance. 



the CAPS has been plotted after averaging multipoles in 11 
intervals. As expected (see Figure 2), all the measured points 
are inside the acceptance interval at 32% significance level. 
We remark, as said in Section 4, that the CSMHW and the 
CCF curves are strongly correlated, whereas the correlation 
along the CAPS values is very small. This indicates that, by 
taking into account all the data and all the correlations in 
the proper way, the three techniques could provide a similar 
detection of the ISW, although differences are expected due 
to the presence of the mask. This will be considered in the 
next Subsections, when constraints in the nature of the dark 
energy are determined. Besides, due to the wavelet proper- 
ties, almost the total signal of the ISW can be concentrated 
in a very specific scale range. 

We have checked that this cross-correlation signal is not 
caused by systematic effects. We focus in the CSMHW de- 
tection, since is the one with the highest significance. In 
particular, the NVSS map has been also cross-correlated 
(independently) with each one of the WMAP receivers (Ql 
and Q2 at 40.7 GHz, VI and V2 at 60.8, Wl, W2, W3 and 
W4 at 93.5 GHz). We found, for all the cases, the same 
CSMHW as for the combined WMAP map. We have also 
tested that the cross-correlations for the maps Q1-Q2, VI- 
V2 and W1-W2+W3-W4 with the NVSS data are compati- 
ble with simulations. For all these noisy maps, the CMB and 
foregrounds contributions are completely removed, showing 
no cross-correlation with the NVSS data. Hence, we con- 
clude that all the WMAP receivers are producing the same 
cross-correlation signal and that it is not due to noisy arti- 
facts. 

We have also checked that the cross-correlation signal 
is not due to foreground emissions. In particular, the ex- 
pected positive correlation between the radio sources pre- 



sented both in the NVSS and the WMAP maps is consid- 
ered. This correlation is expected to show a frequency depen- 
dence, reflecting the emission law of the point source popula- 
tion (e.g. Toffolatti et al. 1998) . The cross-correlation curves 
for the different WMAP frequency channels (Q at 40.7 GHz, 
V at 60.8 GHz and W at 93.5 GHz) are shown in Figure 5. 
We have also considered the cross-correlation curve for the 
map K-2.64Ka, where K and Ka are the lowest frequency 
WMAP channels at 22.8 and 33.0 GHz, respectively. These 
channels are clearly contaminated by synchrotron emission. 
However, an additional CMB map can be generated by sub- 
tracting the Ka map from the K one, multiplying the first 
one by a factor of 2.64. This number corresponds to the 
expected increment of the synchrotron emission from 33 to 
22.8 GHz§. As seen in the left panel of Figure 5, the same 
cross-correlation curve is obtained from the whole WMAP 
frequency range, showing not frequency dependence at all. 
Moreover, the cross-correlation curve for the W-V-Q map 
(free of CMB signal but with a clear foreground contri- 
bution) is perfectly compatible with zero (Figure 5, right 
panel). Hence, we conclude that the cross-correlation is only 
due to the CMB, with a negligible contribution from fore- 
grounds, as expected for the ISW signal. 



5.2 Constraints in Qa 

By comparing the measured detection with the theoreti- 
cal curves for different cosmological models, limits on the 

§ A power law is assumed for the frequency dependence of the 
synchrotron emission: 1 syn 

(u) oc T S yn,{vo){ v I v o) 2 ' 7 i as Pro- 
posed by Bennett et al. (2003b). 
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Figure 6. In the left panel, the likelihoods of the considered statistics (CAPS, CCF and CSMHW) are plotted as a function of Q.a- In 
the right panel we show the cumulative probability of having an Qa greater than a certain value. 



cosmological parameters can be established. In this first 
approach, we are interested in putting constraints on the 
amount of vacuum energy density Oa, assuming a standard 
inflationary universe with w = — 1, without redshift evolu- 
tion. We have explored < Q,a < 0.95. For the bias pa- 
rameter, we have considered the analysis of the NVSS cata- 
logue (using the auto-correlation function) made by Boughn 
& Crittenden (2002, 2004) adopting the Dunlop & Peacock 
(1990) RLF1 model. They found that observations favor bias 
values of b — 1.3 - 1.6. For our analysis (where H = 0.72 
and Q, m = 0.29), b = 1.6 is the best value. In any case, 
we have checked that variations of the bias in the interval 
1.4 < b < 1.8 do not change the results. The Dunlop & 
Peacock (1990) RLF1 model was also used to describe the 
dN/dz distribution since it fits well the NVSS galaxy auto- 
correlation, as it was already shown by Boughn & Crittenden 
(2002, 2004) and Nolta et al. (2004). 

In practice, the observational curves for the CAPS, CCF 
and CSMWH (given by Equations 1, 6 and 8) are compared 
with the theoretical predictions (Equation 2, 7 and 13). For 
each statistic, a generalized \ 2 m defined: 

X 2 (p\b) = A^ ros «,(p|b; Oi)C(a,i, a ] )~ 1 A cross (p\b; a,j) (18) 

where the parameter p = Qa and A cross (a) = 
(cross the ° WN (p\b; a) — cross WN (a)) is the difference be- 
tween the observed CAPS/CCF/CSMHW and the theoret- 
ical prediction at scale a (i.e., the CSMHW R scale or the 
angular scale 9 of the CCF or the multipole £ of the CAPS). 
The matrix C(ai, dj) is the correlation matrix of the WMAP 
and NVSS CAPS/CCF/CSMHW and it is obtained from 
1000 simulations. Assuming that each C(ai,a,j) is Gaussian 
distributed, a likelihood L oc exp(—x 2 /2) can be calculated. 
The likelihood obtained for each statistic is represented in 
Figure 6 (left panel). At laCL, we obtain compatible con- 
straints for the vacuum energy density: Qa = 0.72^0 ^ (for 



CSMHW), Q A = 0.73^-^ (for CCF)^ and Q A = 0.73^ ^ 
(for CAPS). 

On the right panel of Figure 6 we represent (for each 
statistic) the cumulative probability of having an Qa greater 
than a certain value. This quantity provides us an alternative 
method to estimate the significance of the ISW detection. 
These plots shows how the CAPS and the CSMHW are more 
powerful than the CCF for this purpose. There are several 
reasons for understanding this behaviour. On the one hand, 
the CAPS statistic has the important advantage that, apart 
from the mask effect, there are not correlations between mul- 
tipoles, which is always a desirable property to build the 
chi-square, although all the available multipoles contribute 
in a similar way. On the other hand, although the CSMHW 
statistic has correlations among different scales, most of the 
signal is concentrated in a narrow angular range (that, in 
this case, is around several degrees), i.e. the most favourable 
condition to use wavelets. On the contrary, the CCF one 
has correlations among angular scales and, in addition, its 
maximum significance is smaller than the SMHW one. In 
particular, the CAPS provides Qa > 0.1 at a significance of 
« 3.5a, which is the highest detection of the vacuum dark 
energy (based on the ISW) reported up to date. 

The reason of those differences in the significance of the 
detection is due to the presence of a mask. In the case of full- 
sky data sets, a maximum likelihood method applied to the 
three quantities should give very similar results, since each 



^ Notice that the likelihood obtained for the CCF is slightly dif- 
ferent from the one given in Nolta et al. (2004, Figure 4d), al- 
though both provide consistent limits at lcr CL. There are differ- 
ences that can explain it. First, in Nolta et al. (2004) the theoret- 
ical ISW effect is removed from the simulations to estimate the 
covariance matrix. Second, the estimate of the CCF is different 
in Nolta et al. (2004), since each pixel is weighted by taking into 
account the number of unmasked pixels at the original resolution. 
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Figure 7. The logarithm of the marginalized pdf for Qde (upper) and w (middle) are presented. The 2D likelihood is also provided 
(bottom). The errors in the parameter estimation at lc, 2<r and 3c are also plotted. First column shows the limits given by the CAPS, 
the second ones is for the CSMHW and the last one is for the CCF. 



of them can be expressed as a linear combination of any of 
the other two. 



5.3 Constraints in Qde and w 

In the present section, we are interested in putting con- 
straints in the nature of the dark energy. In particular, the 
amount of the dark energy density Qde (and so, the mat- 
ter density, since we impose a flat universe) and the equa- 
tion of state parameter w are studied. As it was done in 
the previous Section, we put limits in these cosmological 
parameters by comparing the measured WMAP and NVSS 
cross-correlation with the theoretical curves obtained from 
different cosmological models. The explored range for the 
dark energy density is < Qde < 0.95, whereas the one for 
the equation of state parameter is —2.0 < w < 0. 

In this situation, a generalized \ 2 can De defined as in 
Equation 18, with parameters p = (Qde,w). As it was said 
in the Subsection 5.2, the elements of the correlation matrix 
are assumed to be Gaussian, hence we can define a likeli- 
hood as L oc exp(— \ 2 /2). A joint pdf of the parameters 
f{Q.DE,w) can be calculated by normalizing the likelihood 
in the whole parameter space. The 2D likelihood is shown 
in Figure 7 (bottom panel). It is straightforward to calcu- 
late the marginalized pdf distributions /(Qde) and /(w). 
These marginalized pdfs are plotted in Figure 7 in the up- 
per (/(Qde)) and in the middle (f(w)) panels. The results 
for the CAPS, CSMHW and CCF are given in columns 1, 2 
and 3, respectively. 

We obtain at la CL: Qde = 0.70t°'.ll and w = 
Qde 



-0.75±°;^ with the CAPS; 

~0-75t° 'to with the CCF ; and DE 
-0.70±gJ with the CSMHW. 



0.70± H and w 
0.651^0 and w 



The limits in the dark energy density are compatible 
with the vacuum energy density calculated in the previous 
section. The constraints in w are less restrictive. The stan- 
dard inflationary concordance model with Qde ~ 0.71 and 
w — — 1 is perfectly compatible within the la CL. As ex- 
plained in the previous Subsection, the differences in the 
results provided by the three statistics arise because of the 
presence of the mask. 



6 CONCLUSIONS 

We have performed the cross-correlation of the WMAP first- 
year data with the radio galaxy distribution traced by the 
NVSS catalogue, in the harmonic, real and wavelet spaces. 
This is the first application of wavelet techniques to study 
the cross-correlations between the CMB and the nearby uni- 
verse. We found a clear cross-correlation at scales in the sky 
around 9 ~ 7° with a significance w 3.3cr, in good agree- 
ment with theoretical predictions. Wavelets provide a bet- 
ter detection level respect to statistics based on spatial and 
harmonics space. This is due to the important advantage 
of wavelets to concentrate, in a particular and narrow scale 
range, almost the total signal of the ISW. On the other hand, 
the CAPS and the CCF provide a similar detection level 
only if the whole information spread along all the multi- 
poles/angles is properly combined and no mask is present. 

Even more, by using the CSMHW, it is not required the 
calculation of different theoretical models, which are neces- 
sary to properly combine the whole information in the CAPS 
and CCF cases. The direct detection of the ISW obtained 
with the CSMHW is independent of the particular values 
for different parameters, as the bias, the point source red- 
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shift distribution, the equation of the estate parameter of the 
dark energy, etc. The ISW detection at w 3.3<r is the highest 
significance level of the ISW reported up to date. By calcu- 
lating the CSMHW for different combinations of the WMAP 
receivers, we have proven that this cross-correlation signal 
is not caused neither by systematic effects nor foreground 
contamination. 

However, the comparison of the measured signal with 
the expected theoretical cross-correlation curves for differ- 
ent cosmological models, is needed to put constraints on the 
cosmological parameters. In particular, we have put limits to 
the amount of vacuum energy density Qa (assuming w = 1), 
obtaining compatible constraints at la for the three differ- 
ent techniques: at laCL, fi A = 0.72±g ^§ (for CSMHW), 
0\ = 0-73±° H (for CCF) and Q A = 0.73^'^ (for CAPS). 
Even more, the CAPS provides Qa > 0.1 at a significance 
of « 3.5(7, which is the highest detection of the vacuum 
dark energy (based on the ISW) reported up to date. We 
have used the bias parameter b = 1.6 already calculated by 
Boughn & Crittenden (2002) for the NVSS data and the 
Dunlop & Peacock (1990) RLF1 model for the dN/dz dis- 
tribution. Our results do not change for bias values in the 
interval 1.4 < b < 1.8. 

We have also tested the ability of the cross-correlation 
to put constraints on the nature on the dark energy by 
studying the amount of the dark energy Q.de and the equa- 
tion of state parameter w. We obtain Q.de = 0.70ljj!20 
and w = -0.75tg;!i (at la CL) with the CAPS, Qde = 
0.70tE!;29 and w = -0-75tg;eo (at la CL) with the CCF and 
n DE = 0.65i°;2o and w = -0-7Qto.il ( at lcr CL ) with the 

CSMHW. Our estimation of the equation of state of the dark 
energy is the first one made through the cross-correlation of 
the CMB and the nearby galaxy density distribution. It pro- 
vides an independent determination from that made by the 
WMAP team (Spergel et al. 2003) using CMB and LSS and 
by other groups (Caldwell & Doran 2004, Melchiorri 2004 
and Corasaniti et al. 2004 and references therein) using, in 
addition to CMB and LSS, data coming from SN-Ia and/or 
BBN. 

We remark that the differences found in the determina- 
tion of the cosmological parameters are due to the presence 
of a mask. In the case of full-sky data sets, a maximum like- 
lihood method applied to the three quantities should give 
very similar results, since each of them can be expressed as 
a linear combination of any of the other two. 

Finally, in this paper we have shown that wavelets are a 
very promising tool for studying the correlation of the CMB 
with LSS data. Other effects like Sunyaev-Zeldovich or the 
contribution from point sources could also be studied having 
an appropriate LSS tracer. 
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